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Abstract. The gauge invariant formulation of Maxwell's equations and the 
electromagnetic duality transformations are given in the light-front (LF) variables. 
The novel formulation of the LF canonical quantization, which is based on the 
kinematic translation generator P + rather then on the Hamiltonian P~, is proposed. 
This canonical quantization is applied for the free electromagnetic fields and for the 
fields generated by electric and magnetic external currents. The covariant form of 
photon propagators, which agrees with Schwinger's source theory, is achieved when the 
direct interaction of external currents is properly chosen. Applying the path integral 
formalism, the equivalent LF Lagrangian density, which depends on two Abelian gauge 
potentials, is proposed. Some remarks on the Dirac strings and LF non local structures 
are presented in the Appendix. 
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1. Introduction 

The LF formulation of Maxwell's electromagnetic theory with both the electric and 
magnetic external currents is interesting at least for two reasons: first as a relativistic 
field theoretical model and second as an introduction to the dual description of QCD. 
The crucial step for such a formulation is the electromagnetic duality transformation, 
which allows one to add magnetic currents into the electromagnetic system with electric 
currents. 

The first attempt in this direction, restricted to a classical theory, done by 
Gambini and Salamo pQ, is partially successful. They have selected two gauge 
invariant independent degrees of freedom and they have formulated the continuous LF 
electromagnetic duality transformation in terms of these degrees of freedom. They 
have shown that this duality transformation is consistent with the usual duality 
transformation for electric and magnetic fields E and B. However when trying 
to incorporate the Cabibbo-Ferrari theory |2j into the LF formulation, they found 
incorrectly 4 gauge invariant independent degrees of freedom. Also Gambini and Salamo 
have proposed effective Lagrangians with the noncanonical current terms. 

Next, some 5-6 years ago, there were two independent attempts to describe the 
discrete electromagnetic duality within the LF description of Maxwell's theory. They 
were initiated by Susskind 3J, who has argued that the LF electromagnetic duality 
transformation can be expressed simply in terms of the transverse potentials as 



This conjecture has been explicitly checked by Brisudova j3j, within the LF canonical 
quantization in the LC gauge A_ = A + = 0, with a final conclusion that only for free 
fields one can define such duality. This follows directly from the starting point, since any 
one gauge potential description of electromagnetic theory is evidently false for a theory 
with electric and magnetic currents. Instead one should use either a two gauge potential 
approach or a LF version of Dirac strings. Also the idea of entanglement of duality 
with gauge transformations [3], jl] is physically misleading - these transformations are 
fundamentally different and should be treated separately. 

The next attempt, by Mukherjee and Bhattacharya [3], is based on a different 
approach with two gauge potentials. This method, introduced long ago by Zwanzinger 
|Hj within the equal-time (ET) formulation, leads to a quite complicated picture, where 
one has to introduce a constant space-like vector n M , n 2 < 0, which breaks the Lorentz 
invariance of starting Lagrangian density. This Lagrangian leads to many constraints, 
which follow from two gauge symmetries. Mukherjee and Bhattacharya, implementing 
the Dirac canonical quantization procedure for systems with constraints 7\, found 
that effectively only one component of each gauge potential is an independent degree 
of freedom. Therefore after redefining these independent modes as two transverse 
components of some effective gauge potential they prove the Susskind conjecture ((TJ) for 
the interacting electromagnetic theory with both electric and magnetic external currents. 
Their analysis of the quantum theory ends with the structure of LF commutators and 



Ai i > SijAj. 
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Hamiltonian. Since their Hamiltonian contains a rotationally noninvariant term, which 
describes instantaneous interaction of electric and magnetic currents, then we may worry 
if a covariant perturbation theory follows from their analysis. 

The aim of our paper is to formulate a clear LF description of Maxwell's theory, 
where only the gauge invariant objects are used. We will start with the free LF Maxwell 
equations and define the LF electromagnetic duality transformation as their symmetry. 
Then we will quantize the system canonically, treating x + as the LF evolution parameter 
and using x = (x~,x 2 ,x 3 ) as coordinates on the LF (hyper-) surface. The Poincare 
generators will be defined in terms of the symmetric energy- momentum tensor Tjf^, 
which is gauge invariant, rather than the canonical energy- momentum tensor T^ n , which 
is a gauge dependent object. We will propose a novel canonical procedure for the LF 
systems, where the generator of translations in x~ direction is used for deriving canonical 
LF Poisson brackets for all independent fields. Since this generator is a kinematic one, it 
will keep its free field form also for an interacting theory, where a Hamiltonian contains 
some interaction part. Our canonical procedure is explicitly duality invariant and we will 
prove the Susskind conjecture (JIJ, but for fields with a different physical interpretation. 

Our paper is organized as following. In section 2 we present a very concise ET 
description of Maxwell's equations with external electric and magnetic sources. In 
section 3 we start with the LF tensor formulation of Maxwell's equations and the 
electromagnetic duality transformation. We also introduce the LF notation for the 
electromagnetic fields. In section 4 we analyze the free field case, when no external 
sources are present. In section 5 we consider the general case of the Maxwell equations, 
when both kinds of external sources are present, paying a special attention to the 
LF duality transformation. In section 6 we switch our analysis to the path integral 
formulation and derive all propagators, which mediate interactions between electric 
and magnetic currents. We also find a local Lagrangian, which is equivalent to the 
proposed Hamiltonian formulation. In conclusions we discuss our results indicating the 
crucial points and suggesting possible further investigations. The LF notation and the 
Green functions are presented in the Appendix A. In the Appendix B we give a short 
presentation of different Dirac's strings for the ET and LF formulations. 



2. Vector notation for Maxwell's equations 

Before starting our novel LF approach, let us briefly review the ET formulation. The 
magnetic charges p m and currents J m are formally added to the electric charges p e and 
currents J e , when one writes Maxwell's equations 

VxB = -^E + J e} V-E = p e , (2a) 

Vx£ = -^B-4 V-B = Pm . (26) 
This set of equations is invariant under the electromagnetic duality transformation 

E i — > B i — > —E, (3a) 
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J e i > J m i > J e , (3^) 

Pe 1 > |»m 1 ► "Pe- (3c) 

This vector notation is not a suitable starting point for the LF formulation, since for 
3-vectors, contrary to 4-vectors, there is no definition of LF components. 

Usually, when there are no magnetic sources (p m = J m = 0), then a pair of 
homogeneous Maxwell's equations (|25jl . may be removed with the help of 4-vectors 
Afj, = (A , A) by the standard definitions of gauge field potentials: 

E = -VA - JU, B = VxA. (4) 

This allows to express the inhomogeneous Maxwell equations (j2oj) in terms of 4-vectors 
A^, which further can be easily converted into the LF notation. This may explain why, 
in almost all papers on the LF quantization of Maxwell's theory, only the gauge field 
potential approach has been used (for reviews see 0, [TP]). 

However if magnetic sources are present, then the formulation in terms of a single 
gauge field potential runs into inconsistencies. The way out is either by means of the 
Dirac string [IT], or by the Wu-Yang potentials to mention only the best known 
solutions. But we would like to stress that the gauge potentials are only quite useful, but 
are not necessary, for a consistent quantization of Maxwell's equations (|2al l26j) within 
the ET approach [T3j, Actually, they couple locally to external currents, but the 
canonical commutation relations can be solely expressed in terms of the electromagnetic 
fields (E,B). Thus we expect that also within the LF approach one can consistently 
quantize Maxwell's equations directly in terms of electromagnetic fields. 



3. LF notation for Maxwell's equations and electromagnetic duality 

As a starting point for the LF formulation of Maxwell's electromagnetism with classical 
external electric sources, we take the tensor formulation: d^F^ v = J M , & vy,p d v F\ p = 0. 
These equations can be easily transformed into the LF coordinates as the inhomogeneous 



LF Maxwell equations: 

<9+£_ = diEi + T, (5 a) 

d + B, = -d_E i -e lj d j B_ + J\ (56) 

= + d.B, + J+ (5c) 

and the LF Bianchi identities: 

d+Bi =d-Ei-diE-., (6a) 

d+B_ = eijdiEj, (66) 

<9-B- = tijdiBj, (6 c) 

where the LF electromagnetic fields are defined as E_ = Ei = F +i , Bi = 



F_i, €ijB_ = Fij. These equations, being manifestly gauge invariant, may be taken 
as the basis for the further canonical procedure, both for the free field case and the 
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interacting field case. 

Next, the tensor form of the electromagnetic duality transformations 

can be rewritten in the LF coordinates as 

E_ i — > B_, B_ i — ► (8a) 

i — > ey-Ej-, A 1 — ► ~ e ijBji (86) 

and hereafter we will refer to them as the LF electromagnetic duality transformation. 
If one doubts that (}8Dj) is correct, then one may introduce another notation: 

E 2 = E 2 , E 3 = —B 2 , B 2 = E 3 , B 3 = B 3 , (9a) 

which allows to express (|8S)l as 

Ei i— > Bi, B t i— > — Ej. (96) 

However we think that the very form of the LF electromagnetic duality is not important 
- what really matters is how these transformations act on Maxwell's equations. One can 
easily convince himself, that in the absence of external sources ( J M = 0) the LF Maxwell 
equations (fo"al ISDl I5c|) and (fflal 1661 KTcf) transform, mutually in pairs, under (|Sal 186)) . Since 
the duality transformation for sources (j36| ETcj) can be directly expressed in terms of the 
4-currents = (p e , J e ) and = (p m , J m ) 

J* , — > k» i — ► -J^, (10a) 

then for the LF components one has 

J* i — >K ± i — ► -J ± , f i — ► tC i — ► -J\ (106) 

At last, we also need a gauge invariant energy- momentum tensor, thus we take the 
symmetric energy-momentum tensor 

T^ m = F\F X » + \g» v F x »F Xp , (11) 

with the LF components: 



T: ym = l -(E*_ + Bl) 1 T+ + m = Bl T^ m = E_B l -e l3 B ] B_ ) (12a) 
T sym = ~ E-Ei- e ijEjB_, T~~ = Et, (126) 
T% m = - {EkBi + E 3 B % ) + 5 tl [E k B k + + l -B 2 _^ . (12c) 

These components of T^ m are invariant under the LF duality transformations (f8~a| l86j) 
and in our further investigations we will always keep our quantization procedure both 
the gauge and duality invariant. 
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4. Gauge invariant canonical quantization for free fields 

In this section we would like to focus our attention on the quantization of electromagnetic 
degrees of freedom, thus we will consider the case of free electromagnetic fields when all 
external sources vanish (J^ = 0). In our LF canonical procedure we choose x + as the 
temporal evolution parameter, thus all classical LF Maxwell equations ([Hal - ItTc j) should 
be classified either as equations of motion or constraints. Usually, when an equation 
contains a term with the temporal partial derivative d + it is classified as an equation of 
motion, otherwise it is a constraint. However here situation is more tricky, since both 
equations (|56| IFTaj) contain terms d + Bi, thus their linear combination give rise to the 
effective constraints 

2d_Ei = d % E_ - e tj djB_, (13) 

and the effective equations of motion 

2d + B i = -d i E„-e ij d j B_. (14) 

Thus we conclude, that for Maxwell's theory there are more LF constraints then ET 
ones. A formal reason for this difference follows from the inhomogeneous equation (j5fejl . 
which being an Euler-Lagrange equation, generally contains the temporal derivative 
term d + Di, where the canonical momentum field is defined as 

A = |f • (15a) 



On contrary, the homogeneous equation (j6o|) , being a Bianchi identity, contains the term 
d + Bi. Generally, these two fields (Di, Bi) form a canonical pair and are independent 
degrees of freedom. However, for Maxwell's theory, with the Lagrangian density 

C = -7 V" = I (E 2 _ - B 2 _) + EtBi, (15b) 



4"^ 2 

one finds that these the canonical variables (Di, Bi) are constrained 

A = B h (15c) 

and two different equations contain the same term d + Bi. Further consequences of this 
observation will be presented elsewhere and here we only stress that the constraint 
(115 cj) is gauge invariant, thus truely a physical phenomenon, which can be intimately 
connected with the Lorentz symmetry of Maxwell's electrodynamics. 

Since there is no equation of motion for E iy then we would like to treat these fields 
as dependent field variables and remove them, from our canonical analysis, by means of 
the constraint equations (fT3^) |. 

All remaining LF electromagnetic fields (£L, Bi) have their equations of motion 

(2d+d- - A ± ) E- = 0, (16a) 
(2d+d- - A ± )B_ = 0, (166) 
2d + Bi = - diE_ - eijdjB^, (16c) 

| This is quite similar to the nondynamical components of fermion field ip- , ip_ , which are removed by 
solving the constraint part of Dirac's equations. 
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and also appear in the constraints§ 

d_E_ + diBi = 0, (17a) 
d_B_ - eijdiBj = 0. (176) 

Further we notice, that only these LF fields appear in the gauge invariant Poincare 
generators of translations 

d 3 x l(E 2 +B 2 _), (18a) 



p- = 


f d 3 x 


T+- = 


sym 


p + = 


1 d 3 x 


rp + + _ 


sym 


p i = 


1 d 3 x 


rp+i 


sym 



2 

d 3 x Bl (186) 

d 3 x (E-Bi — eijBjB-) . (18c) 

Therefore, using the equations of motion (jl6a| I166| 116c)) . one can prove that these 
generators are the constants of motion (do not depend on x + ). 

Since our system contains constraints for canonical field variables, then our next 
steps within the canonical quantization procedure could be taken according to either 
the Dirac method [7j or the Faddeev-Jackiw (FJ) method of the Hamiltonian reduction 
[Bj. However, we prefer to choose a less technical procedure and solve the constraint 
equations (|17a( I176|) by a suitable parameterization of electromagnetic fields: 

E_ = -diAi, B_ = SijdiAj, B t = d„Ai, (19) 

where Ai are the independent gauge invariant fields. This parameterization clearly 
shows, that within the LF formulation, there are only 2 independent dynamical modes 
with the equations of motion 

(2c\<9_ - A ± ) Ai = 0. (20) 

Thus we have a complete agreement between the ET and LF formulations: in both 
approaches we have 2 independent relativistic modes, which are described by 4 canonical 
ET fields and 2 canonical LF fields. 

This observation agrees with pQ, though Gambini and Salamo take E_ and B- 
(using a different notation for them) as the independent electromagnetic modes. Also it 
forms the physical explanation of the additional constraints (J13J) . which appear in the 
LF approach - there is only 1 independent canonical field variable for any relativistic 
independent mode, contrary to the ET approach, where every mode is described by 2 
canonical fields. 

Using our parameterization of fields, we may express the gauge invariant Poincare 
generators as 

i J d 3 x [(diAif + (e^diAj) 2 ] , (21a) 



p- = 


f d 3 x 


T+- 

sym 


p + = 


1 d 3 x 


rp + + 

sym 


pi = 


f d 3 x 


rp+i 

sym 



J d 3 x (d-AiY, (216) 
-J d 3 x (diAjd.Aj) , (21c) 

This is quite analogous to the status of ET electromagnetic fields (E, B). 
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while the LF duality electromagnetic transformations (|8al l8~S)) boil down to 

Ai i-> -eijAj, (22) 

which agrees with Susskind's conjecture ((TJ). Here we notice that our results have the 
same form as those found within the usual gauge potential approach for the the LC 
gauge condition A_ = A + = 0, provided we set Ai = We stress that it is just a 
coincidence, since our variables Ai, being gauge invariant quantities, are true physical 
fields and there is no gauge transformation for them. For free Maxwell's theory, when one 
chooses the LC gauge condition A_ = A + = 0, all nonphysical modes are removed and 
the transverse the potentials Ai describe physical modes. However, for an interacting 
theory, ignoring the difference between Ai and Ai, one may run into serious difficulties 
- just like in jlj. 

Also we would like to indicate another difference between the gauge independent 
field variables Ai and the usual gauge field potential A^. In our case, we use At to solve 
identically two LF constraints ()17al I17fo|) . while the usual vector gauge potentials A^ 
solve identically all Bianchi identities (jHl Eil E3) • 

Since we have effectively reduced the constrained system (E_, B_, Bi) into the 
independent dynamical fields Ai, then the LF canonical procedure may follow directly 
from the Poisson bracket relation 

d+Ai = {Ai,P-} pB , (23) 

when we take the equations of motion (}20|) and the LF Hamiltonian P ()2 1 aj) . However 
here we would like to propose a novel LF canonical procedure, which uses P + instead 
of P~. Thus we argue that one may start with the trivial relation 

d-A = {A,P + } pB , (24) 

and using the expression ()216|) for P + , one may easily infer the canonical LF Poisson 
bracket 

2{d^A l {x + ,x),A J {x + ,y)} pB = -5 tJ 5 3 (x-y). (25) 

Since P + is the kinematic generator, then the relation (|2*4*j) can be effectively used for 
finding the LF canonical brackets also for an interacting theory. 

For the consistency check, one may use the LF canonical brackets (j25j) to calculate 
other Poisson brackets 

&Ai = {Ai,P tl } PB , (26) 
finding (for fi = — ) the proper form of the equations of motion ()2()|) and (for fi = j) the 
trivial identities. 

The canonical quantization procedure means that canonical variables are changed 
into quantum operators and Poisson brackets are transformed into commutators || 

2 [d^Ai(x + , x), Aj(x + , y)] = -%5 3 (x - y). (27) 

| We will denote the quantum field operators by the same symbols as the respective classical fields 
hoping that this will not lead into any misunderstanding. 
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We stress that both the classical brackets and the quantum commutators are invariant 
under the LF duality transformation (}2"2"|) . thus our canonical quantization procedure is 
both the gauge and duality invariant. 

5. Electric and magnetic external currents 

As a next step we would like to consider the case of LF electromagnetic fields 
interacting with electric and magnetic currents. Therefore we take the inhomogeneous 
Maxwell equations with electric external currents (|Hat l5fe| IKcj) and applying the duality 
transformations ()8al \8~R I106|) we generate the equations with magnetic currents 

d + E_ = d t Ei + J-, = tijdiEj + K~, (28a) 

d + Bi = -d_Ei - eijdjB- + J\ d + B t = d^E { - + e i:j K j , (28b) 
«9_£_ = -diBi - J + , d_B_ = eijdiBj - K + . (28c) 

These LF Maxwell's equations with electric and magnetic currents are equivalent to 
()2qJ) and (j2~Sj) and they form the starting point for the analysis of interacting theory. 
There are two consistency conditions for these Maxwell's equations: the electric current 
continuity equation 

d + J + + <9_ J- + dif = 0, (29a) 

and the magnetic current continuity equation 

d+K + + d^K" + dilC = 0. (29b) 

In this paper, we will not analyze the dynamical structure of these electric and magnetic 
currents, thus we will only suppose that these two conservation laws always hold. We 
stress that all equations ()28all28&ll28a|) are inhomogeneous, thus we argue that generally 
no equation can be interpreted as a Bianchi identity and there is no one gauge potential 
approach, which would lead to these equations. 

We will carry our LF quantization procedure following the same steps as in the 
previous section and start with the identification of equations of motion and constraints. 
As before Ei is a nondynamical field variable, which now can be determined from the 



effective constraint 

2d-Ei = diE. - tijdjB- + f - e i6 KK (30) 

The other fields are dynamical with the effective equations of motion 

2d+Bi = -e ij d j B_-d i E_ + J i + e ij K j , (31a) 

(2d + d_ - A ± ) B_ = e^P + diK* + 2d^K~, (31b) 

(2d + d_ - AJ E- = - eijdiRl + dj 1 + 2<9_ J", (31c) 

while the remaining constraint equations are following 

d-E_+diBi = -J+ (32a) 

d-B_ - tifiiBj = - K + . (325) 
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These last equations can be interpreted as the LF electric Gauss law and the LF magnetic 
Gauss law, respectively. Similarly, to the previous section, we wish to find a useful 
parameterization of the LF electromagnetic fields, which would solve identically the 
above constraints. However, now the best choice is no longer evident, since we have to 
allow for some explicit dependence of electromagnetic fields on the external LF charges 
J + and K + . Therefore we may try different possibilities, for example, we may take 
either 

E_ = - diAi - (d-r 1 J + , (33a) 
= SijdiAj - (336) 
Bi =8-Ai, (33c) 

or 

E_ = - diAi, (34a) 
£?_ = eijdiAj, (346) 
B, =d-A i -d i Al 1 J + + e ij d j Al 1 K+. (34c) 

Apparently this arbitrariness is annoying as an ambiguity of our approach, however 
we interpret it as a manifestation of different possible choices of independent modes 
for the interacting system. Every different choice is self-consistent and various choices 
should ultimately lead to the same physical predictions, though they may differ at the 
intermediate stages. Our novel canonical procedure offers some extra hints for a choice 
of parameterization. If we demand that the translation generator P + remains its free 
field forms 

P+ = J d 3 x B? = J d 3 x (d.Ai) 2 , (35) 

then the only acceptable parameterization for B^ is (|33cf . The nonlocal (in x~) 
dependence of electromagnetic fields E- and B- on J + and K + , given by ()33aj) and 
(133 6|) . can be interpreted as constant vector strings. Since these strings appear while 
we solve the Gauss law constraint equations, we may call them the Gauss law strings. 
A more detailed discussion of these strings and their relation to the usual Dirac strings 
is presented in the |Appendix B 



Now it is quite an easy exercise to check that all equations of motion (|31al - I31cjl 
are equivalent to the equation of motion for Ai fields 

(2c\cL - A ± ) Ai = €adj{d-)- x K + + di(d-)- x J + + f + eijK j , (36) 

provided the covariant conservations laws for the electric and magnetic currents (j29aj 
1296)1 are taken into account. Thus we see that in our gauge invariant description 
of Maxwell's theory with electric and magnetic external currents there are still 2 
independent dynamical degrees of freedom, just like in the free field case, though now 
we have also constant string structure (<9_) -1 for the interaction terms. 

The canonical Poisson brackets for Ai fields follow immediately from the trivial 
equation 

8-A = {A u P + } pB) (37) 
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with P + given by ()35|). leading to 

2{d_A l (x + ,x),A J (x + ,y)} pB = -5 l0 5\x-y), (38) 
which is precisely the free field result. Similarly, also P l is a kinematic generator 

al 3 x (diAjd-Aj), (39) 



but P~ , being the LF Hamiltonian, should contain interaction terms. The very form of 
P~ should be consistent with the effective equation of motion (|36j) . which we write as 
the Poisson bracket equation 

2d + d_Ai = 2 {d-Ai, P~} = A ± Ai + e ij d j (d_y 1 K + + ^(d-)" 1 J+ + J 1 + eyiP.(40a) 
Next, due to (|38jl . we infer the functional differential equation for P~ 

- TT = A ± A + e ij d j (d^)- 1 K + + S,^)- 1 J + + J 1 + eyK 3 , (40b) 
which can be solved as 

+ J d 3 x [diAi(d-)- 1 J + - tijdiAjid-)- 1 ^ 

- J d 3 xMJ* + ^ K3 ) + Hcur- (41) 

H cur is a constant of functional integration, which may depend on external currents J M 
and and we will call it the current Hamiltonian. 

We notice that the above canonical procedure is invariant under the LF 
electromagnetic duality transformation, which again has the simple form for the 
independent modes 

A{ i ► tijAj, (42) 

provided the current Hamiltonian in (|41|) is also invariant under the duality 
transformation of external currents (JlOfcj) . Actually, the determination H cur is not a 
trivial task and lies beyond the scope of the canonical quantization procedure. In the 
next section, we will demand a covariant form of the perturbative propagators and this 
will uniquely fix H cur . 

As a concluding remark of this section, we stress that even though the expression 
for the LF Hamiltonian (|41j) is quite similar to the well known LF Hamiltonian for the 
LC gauge A_ = 0, then this similarity is quite misleading - one must be aware that 
there is no way to derive this Hamiltonian within the one gauge potential approach and 
this explains the failure of [3] for the case of interacting Maxwell's theory. 



6. Path integral for electric and magnetic sources 

The canonical structure, which we have found in the previous section, forms a good 
starting point for the canonical quantization procedure for the sector of electromagnetic 
fields. There is yet an undetermined current Hamiltonian H cur , which describes the 
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LF instantaneous interaction of currents and we have to fix it ultimately. We believe 
that the simplest way to find the consistent current Hamiltonian H cur follows from the 
structure of perturbative propagators. This is based on the observation, that within the 
LF formulation, quite frequently one has to supplement a chronological (in x + ) product 
of quantum field operators by some LF instantaneous contribution from a Hamiltonian, 
when one finds a LF perturbative propagator. Since perturbative propagators are 
derived most easily from a path integral definition of a generating functional, therefore 
in this section we will depart from the canonical quantization procedure but rather 
concentrate on the path integral formulation. 

We define the generating functional of all Green functions as the phase-space path- 
integral, which follows naturally from the canonical structure given in the previous 
section^ 

Z[J + , J\K + , K*] = J VAi exp(^i J d A x d+Aid-A^j exp ( -i J dx + P~) . (43) 

Here we stress that the path integral field variables A% are the gauge invariant quantities 
and they are the canonical variables in the unconstrained Hamiltonian phase-space. 
Since all integrals are Gaussian we may immediately perform them and get the result 

Z[J + , J\ K + , K l \ = exp l - J d i xd i yK^(x)g flu (x - y)K v (y) 
x exp l - J <Pxd t yJ' i (x)g ia ,(x - y)r(y) 



x expi j d 4 xd 4 yJ t "(x)V flu (x-y)K I/ (y) 

x exp ( -% / dx + H cur ), (44) 



with the propagators 

Gpv = (-<//«/ + O/A + n„0 M )(n • d) -1 ) D F - n^n^n ■ d)~ l (n ■ <9)~\ (45a) 
V, v = t, vXp n x dP{n-d)- 1 D Fl (456) 

where n M is the null vector (n+ = l,n = 0), thus (n ■ d)~ x = (cL) -1 . The instantaneous 
part in (|45aJ) can be removed if we choose the current Hamiltonian as 

H cur = \jd 3 x [(9_)-V + ] 2 + \ J d 3 x [(d^K+f . (46) 

This leads to the effective propagators for the electric-electric and magnetic-magnetic 
sectors 

2V = {~9nu + (n^ + n„d^)(n ■ d)' 1 ) D F , (47) 

which have the form of the Abelian gauge field propagator for the LC gauge condition: 
n^A^ = A + = A_. Our propagators depend on the null vector n M , which enters into the 
theory due to the LF quantization procedure. However, when the external currents are 
conserved, then the dependence on n M disappears for the propagator but remains for 

% We will omit the normalization constants for all path-integrals which will appear in this section 
keeping in mind the normalization condition Z[0] = 1. 
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the propagator T>^ u . This means that the Lorentz invariance is broken by the quantum 
interaction of electric and magnetic currents. This is a real challenge to prove that for 
the physical observables (like cross sections etc.) one may restore the Lorentz symmetry. 
Within the ET formulation, one may prove |18j . [T9*] . [2U] . that the dependence on a 
space like vector n M finally disappears. We hope that a similar phenomenon happens 
also for (J37J), though here we have a dependence on a null vector n 2 = 0. 

The Fourier transform of these propagators contains the Principle Value (PV) 

prescription for the noncovariant pole PV- — . It is known, within the ET formulation, 

k_ 

that this prescription is not consistent for a nonAbelian gauge field theory |17j . In 
our gauge invariant procedure, we see that the PV prescription is the only possibility 
for solving the LF constraint equations in terms of a real valued distribution, at fixed 
x + . Thus the status of the PV prescription is consistent here within the LF canonical 
formulation of Maxwell's theory. 

The form of our perturbative propagators strongly suggests, that there should be 
some gauge field model, which effectively produces the same path-integral as we have 
found here. Since there are non local interaction terms in the current Hamiltonian ([46)1 . 
thus we may add two auxiliary path integral variables A + and C + into the phase space 
path- integral (|43|). This allows us to write the equivalent expression for the generating 
functional 

Z[J + ,J i ,K + ,K i ] = J VAi VA + VC+ x 

x exp i J d 4 x [Ci oca i + C+K+ + A + J+ + A(f + eijK 1 )] , (48a) 
where the local Lagrangian density C-i oca i is 

Clem = \ (d.A + ) 2 + (d+A, - d t A + ) d.A t + i (<9_C+) 2 + e lJ d l A J d.C + . (486) 

In these expressions we have changed our gauge invariant path- integral variables Ai into 
Ai, since here they are just dummy variables of path integrations. 

In next step we treat this local Lagrangian density Ci oca [ as the case of the double 
LC gauge condition A_ = C_ = imposed on the gauge invariant Lagrangian density 

Cinv = \ (d + A_ - d^A + f + (d+A, - d,A + ) (9_A t - d t A_) 

- ~ {eifrAj) 2 + ~ (d + C_ - d_C + - t i3 d % A 3 f , (49) 
where we have two independent Abelian gauge transformations: 

Ap — > A^ + d^Qe, (50a) 
C ± ^C ± + d ± Q m . (506) 

The first three terms in (|49J) form a standard LF Abelian Lagrangian with the gauge 
potentials A^ while last term describes the contribution of the dual potentials C±. We 
stress that without this last term one cannot consistently quantize Maxwell's theory 
with both electric and magnetic currents. 
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Since in previous sections, we kept our canonical procedure explicitly duality 
invariant, thus we would like to check whether ()49|) has some duality symmetry. We find 
that, this Lagrangian density is invariant (up to a total derivative) under the following 
transformation 

Ai i— > -UjAj + (d-)- 1 (diC- + tijdjA.) , (51a) 
A± i — ► C± i — > -A±. (516) 

Since for the double LC gauge condition A_ = C_ = 0, these transformations reduce 
to our former gauge invariant duality transformation (|42J1. thus we argue that they 
are the generalized LF electromagnetic duality transformations for the gauge and dual 
potentials. 

With the transformations (j51«j) and ([51 5)1 . supplemented with the transformation 
of currents fjlOajl . we can check that the duality symmetry (up to a total derivative) is 
possessed by the interaction Lagrangian density 

C mt = A_ (j- - e^cL)- 1 ^') + C_ (k~ + {d^d.K 1 ) 

+ A + J + + C + K + + Ai [f + €ijK j ) . (52) 

Apparently this interaction Lagrangian is not satisfactory, since the gauge and dual 
potentials couple nonlocally to external currents. We could introduce further auxiliary 
variables to keep all potential-currents couplings local but at the price of introducing 
the Lagrange multiplier fields - this would be another manifestation of the canonical 
constraint ()15cj) . However we have decided to not proceed in this direction here and to 
discuss it elsewhere. 

Though the interaction Lagrangian (|5*2*j) looks strange, it behaves properly under 
the gauge transformations. From the transformation (j50aj) one gets the electric current 
conservations law = 0, while from (|50 fej) one gets the magnetic current conservation 
law d^K^ = 0. This convinces us that we may treat the gauge invariant Lagrangian 
(I49J) with the interaction Lagrangian (|52j) as a good starting point for the canonical 
quantization procedure with different gauge conditions [2*T] . 

7. Conclusions and further prospects 

In this paper, we show that the electromagnetic duality can be consistently defined for 
the LF Maxwell theory with the classical external electric and magnetic currents. For 
practical reasons we introduce the independent physical fields Ai which allow to satisfy 
identically the LF electric and magnetic Gauss law equations. Remembering that fields 
Ai are not the transverse components of some gauge potential, we are satisfied that for 
these physical modes the LF duality transformation looks like the Susskind conjecture 

We propose a novel canonical LF procedure, which is based on the longitudinal 
translation operator P + , which, being a kinematic generator, has the same form both for 
free and interacting theories. Usually the canonical procedure starts with a Lagrangian, 
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but in the case of Maxwell's theory with electric and magnetic currents, the starting 
point is Maxwell's equations. In our novel canonical procedure, we can safely take P + 
from the free field theory, while the LF Hamiltonian P~ can be found as a solution of 
the functional differential equation. 

We also use the path-integral formulation for achieving two different goals. First, we 
find the perturbative propagators with the instantaneous (in x + ) terms, which then can 
be cancelled by the contribution from the LF Hamiltonian. The effective propagators 
(|45 3 l4*7jl have the form known from Schwinger's source theory ^H] for an arbitrary 
constant vector n M . Here we present their LF canonical derivation for the case of a null 
vector n 2 = 0. Second, we prove that our gauge invariant formalism is equivalent to 
some gauge theory with two potentials. We show that the LC gauge condition can be 
chosen for both the gauge and dual potentials. Our two gauge potential Lagrangians (J49|) 
and (|H2*|) are far simpler than the one proposed long ago by Zwanziger |Hj and recently 
used for the LF quantization by Mukherjee and Bhattacharya j^j- Here we would like 
to point out that our result disagrees with jH], where the Hamiltonian contains an 
additional instantaneous term (for the interaction of electric and magnetic currents), 
which explicitly breaks the rotational symmetry. In our case the rotational symmetry 
is broken solely by a choice of the LF surface. 

The two gauge potential Lagrangian densities ()49|) and (|52*j) can form a starting 
point for the dual formulation of Abelian theory within the LF approach. When 
further such formulation is generalized to a nonAbelian theory, then a very encouraging 
possibility arises - the LF version of the dual superconductor models of color confinement 
(for introduction and earlier references see [22])- 

When one performs the perturbative calculations, with currents generated by 
some charged matter fields, then inevitably the perturbative propagators (j45M l4*7j) 
will produce momentum integrals with the UV divergences. Therefore one will have 
to regularize the model without spoiling the gauge and duality invariance. Since the 
antisymmetrical symbol appears explicitly in many expressions, then the dimensional 
regularization seems to be impractical here and one should look for a kind of the Pauli- 
Villars regularization. Possibly, this can give rise to the regularized Lagrangian density 
with higher derivative terms (see [23 J. 

When the charged matter is treated as a part of the quantum dynamical system, 
then one may check the second Susskind conjecture j^j, that in the LF approach there is 
no need for any other non localities, like Dirac strings, beside the usual (cL)^ 1 integral 
operator. Such investigations should also answer another question how the Dirac- 
Schwinger quantization condition for electric and magnetic charges [H], [23] appears 
within the LF formalism. We hope to give definite answers to these problems in our 
future publication. 
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Appendix A. The LF notation 

We use the natural units c — H — 1. Our LF notation starts with the definitions of null 
components for the coordinates x ± = (x° ± x 1 )/y / 2, while the transverse components 
3 ). The similar definitions are taken for any 4- vectors. The LF surface 
coordinates are denoted ). The partial derivatives are taken with respect 

to contravariant coordinates, thus we have d + = <9/<9x + ,<9„ = d/dx~,di = d/dx i . The 

metric tensor has non vanishing components g + - = g (. = l,gij = —b~ij- The scalar 

product of 4-vectors is a ■ b = a + 6_ + a_6_ — aj&j, while for the LF surface components 
we have a ■ b = a~b- — ajfrj. There are two natural antisymmetric tensors = — €ji 
with e 23 = 1, and e + _jj = e^. 

The inverse differential operator (cL) -1 is taken as the distribution 

(d-)-\x) = ^sgn^-)^^), (A.l) 

which means that we use the PV prescription in its Fourier integral representation 

r rl 3 k - - 1 

(d^-\x) = -if A^^*PV- (A.2) 
Another integral operator djAj 1 is given as 

"1 i 

diA^tx) = —^S(x-), x\ = xV, (A.3) 
2tt x ± 

while the covariant propagator function is defined as 

Appendix B. Dirac's strings within the ET and LF formulations 

Dirac's approach with strings [TT], is basically a gauge potential formulation, where 
the modifications in Bianchi identities (2b) are compensated by the redefinition of the 
electromagnetic field tensor 

= d^A v - d u A^ + ^ uXp G Xp , (B.l) 

provided the string tensor G xp satisfies the equation 

d\G xp = K p . (B.2) 
A straight line string is given by 

G xp = (n • d)- 1 (n x K p - n p K x ) (B.3) 

where n p is a given fixed vector. In the ET approach one chooses this fixed vector to 
be space-like n p = (0, n), thus the integral operator (n ■ d)^ 1 = — (n ■ V) does not 
include temporal evolution. The orientation of n is arbitrary and no physical result 
should depend on it. 
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This procedure can be also applied in the LF approach and one finds that 
there are 2 different possibilities - one can take either (n~ = 1, n + = n % = 0) or 
(n + = n~ = 0, n l ^ 0) - which lead to no temporal evolution in the integral operator 
(n ■ d)^ 1 . In the first case one has a fixed null vector n 2 = 0, while in the second case 
one again has a space-like vector n 2 < 0. 

However we stress that equation (jB.2|) is not duality invariant since it contains 
only the magnetic sources K p . Since our present paper is devoted to the gauge and 
duality invariant formulation we will not discuss such Dirac strings any longer here but 
postpone a more detailed analysis to a future publication. 

Within the Hamiltonian formulation, one can adopt Dirac's idea (|B.1J) for 
transforming inhomogeneous constraint equations into homogeneous equations. One 
may define the electromagnetic field tensor as 

Fpu = + n»U - n v fft + e^ Xp n x g p , (B.4) 

where the components of T^ v satisfy the homogeneous constraint equations, while / M 
and g^ are some functions which may depend on electric and magnetic charges. This 
modification is duality invariant, provided one introduces the duality transformation for 
/" and 

g * _> f» _> - g », (B.5) 

The fixed vector is quite arbitrary and here we will discuss only the simplest choices. 
In the ET approach one can take (n° = 1, n — 0) which leads to 

E = £ + f, V-/ = J°, (B.6) 
B = B - g, V-g = -K°. (B.7) 

In the LF approach one can take (n~ = 1, n + = n % = 0) which leads to 

E_=£_- f + , d_f + = J + , (B.8) 
B_=B-+g + , d-g + = -K + , (B.9) 
Ei =£i, Bi = Bi. (B.10) 

These two choices of a constant vector n p are equivalent, since they indicate the 
respective temporal evolution parameters: in the ET case we have = x°, while 

in the LF case we have n^x^ = x + . Evidently equations (jB.8HB.10|) are the solutions 
(IHHaMHHcj) proposed in the main text, while equations (jB.6HB.7j) lead to a Coulomb-like 
solution for electric and magnetic monopoles. 
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